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Abstract 

We study kinetics of diffusion-limited catalytically-activated A -\- B ^ B reactions taking place 
in three dimensional systems, in which an annihilation of diffusive A particles by diffusive traps B 
may happen only if the encounter of an A with any of the Bs happens within a special catalytic 
subvolumen, these subvolumens being immobile and uniformly distributed within the reaction bath. 
Suitably extending the classical approach of Wilemski and Fixman (G. Wilemski and M. Fixman, 
J. Chem. Phys. 58:4009, 1973) to such three-molecular diffusion-limited reactions, we calculate 
analytically an effective reaction constant and show that it comprises several terms associated with 
the residence and joint residence times of Brownian paths in finite domains. The effective reaction 
constant exhibits a non-trivial dependence on the reaction radii, the mean density of catalytic 
subvolumens and particles' diffusion coefficients. Finally, we discuss the fluctuation-induced kinetic 
behavior in such systems. 
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I. INTRODUCTION 



Catalytically-activated reactions invo 
cesses in physics, chemistry and biology 
can be written as 



ving diffusive species underly many different pro- 
l|, 0,0, J, 0]. In a general notation, such reactions 



A + B + C ^ P + C, (1) 

where A and B designate two different types of mobile reactive species, C denotes a catalytic 
subvolumen, while P stands for the reaction product. Catalytic subvolumens C may form 
some patterns, be uniformly or regularly spread along a given structure, (e.g., a polymer 
or polymers in solution), scattered uniformly within or on the boundary (surface) of the 
reaction bath. 

The reaction scheme in Eq. signifies that the bimolecular reaction between the A and B 
molecules may take place, at some finite elementary reaction rate k, if and only if a diffusive 
encounter of an A and a B happens within any catalytic subvolumen C. In some cases, the 
B particles may be unaltered by the reaction or their concentration may substantially exceed 
that of the A species, in which situation the reaction in Eq.(0) can be viewed as bi-catalytic: 
that is, one deals with a simplified reaction scheme A ^ P which requires the presence of 
two different catalytic subvolumens - B and C. This is most often the case in biology, as 
exemplified, for instance, by the transcription of genes induced by the simultaneous presence 
of several transcription factors on the promoter sequence 0]. 

Most of analytical descriptions of the catalytically-activated reactions have focussed so 
far on the particular question how reactions are promoted by specific catalytic subvolumens, 
which was believed to be the most crucial aspect of the problem 0]. Within this line of 
-hought, the kinetic behavior has been determined using standard formal-kinetic approaches 
3, 0. On the other hand, a few available analytical studies of the catalytically-activated 
reactions limited both by diffusion of species and by the condition of the simultaneous 
encounters within the catalytic subvolumens have revealed a non-trivial kinetic behavior 
in low dimensional systems, and showed that although in three dimensions kinetics follows 
standard temporal behavior, the effective reaction rates are strongly dependent on particles' 
diffusion coefficients, subvolumens' radii and their concentration 3|,y,|5|. These findings are, 
of course, in an apparent contradiction with the predictions of the formal-kinetic approach. 

As a matter of fact, it has been already realized that for reactions taking place in non- 
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catalytic systems, in many instances, the kinetic behavior can not be adequately described 
in terms of the formal-kinetics approach. Indeed, it has been known for a long time that 
diffusion of reactive species limits the reactive process and leads to unusual kinetics in 
low dimensional systems (see, e.g., Ref.jsl and references therein). Moreover, it has been 
discovered that in many reactive systems spatial fluctuations in particles concentrations 
dominate the long-time evolution and entail anomalous, fluctuation-induced behavior. In 
particula... a p.o^ounced devotion from the conventional descriptions Q has been predicted 
for the irreversible, diffusion-controlled recombination reaction A + B ^ in case when 
initially the particles of the A and B species are all distributed at random, independently 
of each other and with strictly equal mean densities ri^(O) = nsiO) = Uq. It has been first 
shown that here at long times the mean particle densities follow 

n(t)~nf (2) 

where d is the space dimensionality and D - the sum of particles' diffusion coefficients, 
D = Da+Db- This law, which was rigorously proven in Refs. should be contrasted to 



the conventionally expected Smoluchowski-type form n{t) ~ l/(j)^^\t) j^, where, as t ^ oo 



' ^^/Dtj-n, d = 1, 



^"Xt) = I drKsir) ~ < 

^ AnDRt, {k = oo) d = 3, 



ln(4Dt/i?2)' 



Ks{t) being the (i-dimensional Smoluchowski-type constant, defined as the flux of diffusive 
particles through the surface of an immobile sphere of radius R - the reaction radius. Note 
that both decay laws contradict to the text-book formal-kinetic description based on the "law 
of mass action", which predicts that regardless of the spatial dimensionality n{t) decays as 
n{t) ~ l/t 0. 

Therefore, according to the Smoluchowski approach, in diffusion-controlled recombination 
reaction A + B ^ Q diffusion slows down the decay in low dimensional systems and entails 
the renormalization of the reaction rates in three-dimensions. On the other hand, in the 
particular case when initially the particles of the A and B species are all distributed at 
random, independently of each other and with strictly equal mean densities 71^(0) =71^(0) = 
riQ, fluctuations in spatial distributions of the reactive species appear as the most important 
rate controlling factor in spatial dimension d < A and dominate the long-time kinetics. 



For trapping reactions A + B ^ B in completely catalytic systems two opposite limiting 
situations were most thoroughly studied. Namely, the case when the A particles diffuse 
while the traps B are static, and the situation in which the A particles are immobile while 
-Bs diffuse - the so-called target annihilation problem (TAP). The case of static traps has 
attracted most of interest prompted by, in part, an early observation [3l that the long-time 
survival probability -Pyi(^) of A particles diffusing in the presence of randomly placed (with 
mean density [B]) traps exhibits highly non-trivial, fluctuation-induced behavior of the form 



lnPA(i)~-[fi] 



2/(d+2) 



(DAt) 



d/{d+2) 



t 



OO, 



(3) 



which stems from the randomness of B distribution in space and namely, from the presence 
of large spatial regions devoid of traps where the A particles survive anomalously long 
times. This fluctuation-induced decay law is intrinsically relevant to the so-called Lifschitz 
singularities near the edge of the band in the density of states of a particle in quantum 



Lorentz j first noticed in 1131. Later works (see, e.g., Refs 



.Q,liilQ,lii 



li 
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have also pointed out relevance of the issue to the problems of 
percolation, self-avoiding random walks or self-attracting polymers, as well as anomalous 
behavior of the ground-state energy of the Witten's toy Hamiltonian of supersymmetric 
quantum mechanics [l4|. 

Survival probability Ptargetit) of an immobile target A of radius R in presence of point-like 
diffusive traps B - the target annihilation problem (TAP), allows for an exact solution in 



any spatial dimension 



24 



Ptargetit) =exp(^-[B]^^^\t)^, 



(4) 



where 4'R\t) has been defined in Eq.® in which one has to set Da = 0. Extensions to 
systems with hard-core interactions between traps [^^l or fluctuating chemical activities of 
traps 



28| have been also provided. 
Now, the general and physically most important case of trapping reactions when both 
A and -Bs diffuse with diffusion coefficients Da and Db was not solved exactly up to the 
present time. It has been proven that in this case the A particle survival probability 
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obeys 

f t'/\ d = 1, 



In PA{t) = -\d{DA,DB) X <^ 



^ d = 2, (5) 



ln(t) ' 
t, d = 3 



which defines its time- dependence exactly. 

On the other hand, the factor Xd{DA, Db) remained for a long time an unknown function 
of the particles' diffusion coefficients and spatial dimension. Only very recently some rigorous 
arg uments have been presented showing that ioi A + B ^ B reactions in low dimensions j^, 
isfll. I31I. . the long time decay of A particles concentration is exactly as in the TAP problem, 
Eq.(j3]), and thus is asymptotically independent of the A particles diffusion coefficient even 
in the case when both species diffuse. 

In three dimensions, however, the precise form of X^iIDa, Db) is still undetermined and 
still very little is known about it. It has been shown that \d{DA, Db) is less than the 
rate constant calculated within the Smoluchowski approach and moreover, it has been 
realized that in case when Da and Db are sufficiently small, Xd(DA, Db) may be bounded 
by a non-analytic function of particles' diffusion coefficients [33|. A perturbation theory 
approach for calculation of \(i{DaiDb) has been proposed in Ref. 2^ and the corrections 
to the predictions of the Smoluchowski approach have been evaluated. It was also shown 
that \d{DAi Db) can not be represented as the function of the combination D = Da + Db 
only, since the diffusion-reaction equations are not separable. Therefore, even in completely 
catalytic systems the evaluation of \d{DA, Db) represents a fairly complicated many-body 
problem. 

In this paper we study the kinetics of the catalytically-activated diffusion-limited reactions 
in Eq.(^ in the special case when B particles remain unaltered by reactions, i.e. the case of 
diffusion-limited catalytically-activated trapping reactions, description of which poses such 
serious technical difficulties even in the non-catalytic systems (see the discussion above). In 
order to obtain an effective reaction rate for such bi-catalytic reactions taking place in a 
homogeneous three-dimensional medium, we first develop an analytical approach, inspired 
by the work of Wilemski and Fixman jSJ], which allows one to estimate the reaction rate 
for non-catalytic bimolecular reactions. Here we extend this Wilemski-Fixman approach 
(WFA) to catalytically-activated trapping reactions. 

We assume that the catalytic subvolumens C are immobile and are spread uniformly in the 



reaction bath. On the other hand, A and B particles are assumed to perform unconstrained 
diffusion, and react only when they are simultaneously present in a spherical domain of 
radius R centered around each catalytic subvolumen C. As well, we suppose that there is no 
other interaction between the particles except for the reaction, which enables us to describe 
the A particle dynamics in terms of a Fokker-Planck-type equation with a sink term which 
mimics the presence of the traps B and of the catalytic subvolumens C. In order to obtain 
a closed equation, we follow then the well known Wilemski-Fixman approach (WFA) js^. 
This approximation relies on the time and space separation of the joint probability density; 
that is, the probability density is assumed to be a product of the equilibrium density and 
of a certain time- dependent function, supposing that initially the system is at equilibrium. 
We hasten to remark, that, although the validity domain of this approximation is still not 
really known, many researchers have shown that the WFA describes quite correctly reaction 



kinetics in several general situations 
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36|. For example, Doi j35| showed that the WFA can 



be used for the purely diffusion-limited case {k = oo), contrarily to the intuitive expectation 
that the WFA is appropriate for systems with a weak chemical reaction rate k. We will show 
also in what follows that the effective reaction constant obtained within such an approach 
reduces to the well-established results in several limiting cases. 

Finally, we will present an estimate of the impact of the fluctuation effects on the long- 
time kinetics of the diffusion-controlled catalytically-activated reactions. 

The paper is outlined as follows: In section II, we formulate the model and write down 
basic equations. Section III is devoted to the solution of the evolution equations within the 
framework of the suitably extended Wilemski and Fixman approach. Here, we determine 
an effective reaction rate for the bi-catalytic reaction and show that it can be is expressed 
through different functionals of Brownian motion, known as residence times and joint res- 
idence times of Brownian particles in some specified domains. Further on, in Section IV 
we calculate the residence times involved and obtain an explicit expression for the effec- 
tive rate constant describing diffusion-limited catalytically-activated trapping reactions. As 
well, we discuss its asymptotical behavior in several limiting situations. Next, in Section V 
we present some estimates of the long-time fluctuation-induced behavior in such systems. 
Finally, in Section VI we conclude with a brief summary of our results and discussion. 
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II. MODEL AND BASIC EQUATIONS 



Consider a three-dimensional reaction bath of volume V comprising a single A particle, m 
traps B, and q immobile catalytic subvolumens C (see Fig^). These catalytic subvolumens 
are uniformly distributed in the reaction volume with mean concentration [C]. The A and 
B particles diffuse freely with diffusion coefficients Da and Db, respectively. For simplicity, 
we assume that A and Bs are point-like particles of zero radius such that excluded-volume 
interactions between them can be safely neglected. Now, the reaction between particles - 
an annihilation of the A particle by any of the Bs - takes place with a given probability 
defining the elementary reaction constant k when both species appear simultaneously within 
a spherical region of radius R (which can be also thought off as the reaction radius) centered 
around any catalytic subvolumen C. Such a "reactive" situation is depicted in Fig|21 

We now proceed by suitably extending the celebrated approach devised originally by 



Wilemski and Fixman 



3J] for non-catalytic trapping A + B ^ B reactions. The basic 



idea behind this approach is that the presence of traps B can be effectively modelled by 
introducing a sink function S into the diffusion equation describing dynamics of the A and 
B particles. Then, the sink term determines the efficiency of the reaction as a function of 
the instantaneous separation distance between the A particle and the traps B. 

In the simplest formulation, this sink term can be represented as the Heaviside function, 
which implies that the reaction takes place with an elementary reaction constant k as soon 
as the A particle appears in the vicinity of any of the traps B. Adapting this line of thought, 
we describe the A and B particles dynamics in terms of the following multivarient diffusion 
equation: 



i 

where \l/({r},t) is time-dependent m + 1 particles probability density function, {r} = 
{r^,rBj, ...,rBm} defines the positions of the A particle and all m traps B, while the sink 
function S is represented as a set of Heaviside functions H{x), {H{x) = 1 if x > and 
H{x) = for X < 0), centered around q catalytic subvolumens: 
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y/2 

i=l k=l 
^ m q 

i=l fc=l 

4 

in which l^' denotes the volume of a "reactive" domain, V' = -vri?^, (see Fig|2)). 

3 

Consequently, the desired probability P{t) that the A particle survives up to time t obeys: 



P(t) = I d{r}^i{r},t), (8) 
and can be readily evaluated once \E'({r},t) is known. 

III. SOLUTION OF THE EVOLUTION EQUATION. 

Equation ^ can be cast into an equivalent form by using the Green function 
G{{r},t] {r°},t°}of equation (jH)) without the sink term, the latter being considered as an 



inhomogeneity 



I . In doing so, we find that the formal solution of equation reads: 



^{{r},t) = I d{r'}vl/({r'},0)G({r},t;{r'},0) 
-k df j d{r°}G'({r}, t; {r^}, tO)5({rO})^({rO}, (9) 
Note that G'({r},t; {r°},t'') of equation © without the sink term factorizes: 



G({r}, t- {r°}, t°) = Ga{va. t; r° , t°) J] GsirB^.t; r%^,t'). (10) 

n=l 

Now, supposing that at t = the traps B were uniformly distributed in the reaction bath, i.e. 
that \l/({r'},0) = ^(0), we find that the first term on the right-hand-side of Eq.Q reduces 
to ^(0). Following the reasonings of Wilemski and Fixman, we assume further on that ^1/(0) 
is the equilibrium density \l/eg = This corresponds to the physical situation when 

the system is first brought to equilibrium and the reaction is triggered then at time t = 0. 
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Next, multiplying both sides of Eq.ijHJ by S'({r}) and integrating it over all spatial variables 
{r}, we obtain: 



t 





v{t) = v,,-kj df J d{r} J d{r'}S{{r}) 

xG({r}, t; {r°}, t°)5({r°})M/({r°}, t°) (11) 



where, by definition, v{t) = J d{r}S{{r})'i/{{r},t) and, equivalently, Veq = 
f d{r}S{{r})'^eq- Using next the definition of ^dr}) given in equation ((Tj), we obtain 
for Vp 



-'eq- 



11""' 



lit ^ n 



^1 ya ym+i 

i=l k=l 

1 1 ft^"-^'^-^- 



YI2 ym+1 y2 
i=l k=l 

Further on, integrating equation (jH)) over all {r}, and using the definition of the survival 
probability (jHl), we find that: 

^^-M.). (1.) 

Consequently, the function v{t) determines the rate of the time evolution of the A particle 
survival probability. 

In order to obtain a closed equation for v(t), we follow again Wilemski and Fixman 
method, assuming that \E'({r},t) can be split into the product of a time-dependent function 
and the equilibrium density corresponding to the situation without reaction, knowing that 
initially \E'({r},0) = \&eg- This approximation is valid a priori for a small enough value of 
the reaction rate k, such that the probability density can be thought of being close to the 
equilibrium's one at any time. As a matter of fact, this approximation still holds in much 
more general situations, for example when the reaction is diffusion-limited {k = oo), as it 



was shown by Doi j^. A detailed discussion of this approximation can be found in Refs. 

Applying it to our case of the catalytically-activated diffusion- limited trapping 
reactions, we have: 
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vl/({r},t)^vl/^,i.(t). (13) 

Noticing next that 

j d{r}S^ = u{t) j d{r}S^eq = y{t)veq = v{t), 
we get uit) = v(t)/veq and hence, the approximated probabihty density reads: 



* = (14) 

Veq 

The latter equation, within the framework of the Wilemski-Fixman approximation, yields 
the following result for v{t): 



v{t) = - f dt\{f)I, (15) 

'^eq Jo 

where the integral / is given explicitly by: 

^ = Vr.Y. f ^{r} / 4r°}i^A,c.i^B.c.i^A0,c,i/BO,c,G({r},t;{r°},t°). (16) 

i,j,k,i 

In order to obtain an explicit expression for this integral, and thus to get an access to 
the kinetic behavior of the survival probability, we consider it in more detail below. We 
note that I can be split into four parts, when specifying the following different events: the 
A particle, initially present in the k-th subvolumen C together with the z-th particle B, will 
further encounter either the same B particle {i = j) or another B particle {i 7^ j) in either 
the same subvolumen C{k = I) or some other subvolumen C {k ^ I). 

First, let us consider the integral lijki for i ^ j. Integrating over r^^ for m 7^ i and over 
r^^ for V ii we obtain: 

Ym-2 r 

hjkl;i=ij 

- -^JT / drAdr°AdrBA%HA,c,HB,,c^HAo,c,HBO^CiGAirA,t;r%t°). (17) 
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Next, performing the integration of the Heaviside functions H^o q^ and HBi,Ck over rfr^- and 



dr^., we have: 



hjkl;ijtj 

= / dvAdr^AHA,c,HAO,c,GA{rA,t;r%f). (18) 

Let us now consider the integral lijki for i = j. Integrating over the variables which do 
not appear as arguments of the Heaviside functions, we can write that: 



^m— 1 

^iikl ~ 



dvAdr^AdrBi dr% HA,Ck Hb^ ,c, Hao ,Ci Hsf ,c, 



xGAivA, t; r% t°)GB(rB„ t; r%,t'). (19) 

Note that when k = I, HA,CkHB,,CkHA»,CiHBf,Ci reduces to HA,CkHBi,Ck- 

Finally, using the latter decomposition of the integral J, and once again splitting it with 
respect to cases k = I and k I, summing over i,j, k, I we obtain: 



/ = m{m - l)qlijkk;i^j + mqlukk + m{m - 1) 

xg(g - + mq{q - l)Iiiki-k^i- (20) 

Now we have to average the integral / over the positions of the catalytic subvolumens C. 
Performing such averaging and denoting it by angle brackets with the subscript r,^, we have 
the following four terms: 

Ym-2 r 

(1) {Iijkk;i^j)vc^ = y,2 j drAdY\HAHAoGA{rA,t]v\,t^) 

Ym.-2 yyl2 ym-2 

(2) {Iijki,i^j,kj^l)rc^,rci ~ ~~yi2 y2~ ~ y 

(3) {Iiikk)rc^ = -yJT j dvAdv\dvB4v% 

ym—l 



rm—l 



Kl(t-t°) 



ym-1 ^ r 

(4) {I^^kl■,k^l)rc,,rc, = -y^y J dudr Adr^Ad^ bA% 

xHaHb^HaoHbo xG^(rA, t; + u, t°)GB(rB^, t; r% + u, t°) 



ym-l 

(21) 
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where x(^); '^i(^) and K2{t) are some functionals of time, which will be made explicit in the 
next section. 

Using the latter results, we may express the v{t) function in terms of the functionals 
Xit - t°), Ki{t - t°) and K2(t - t°). This gives 



v{t) = v^.-k ^'^ ^^^^ j\t%{e)-k'^ j\Mt')x{t-t') 



^ ^ dMt')K,{t-t^)-k\-^ I dt\{f)K2{t-f). (22) 



y.j/Mt')n.{t-t^)-k'LA 
Now, defining the Laplace transformation over the time variable of a function f{t) as 





/(s) = dtf(t)e and performing the Laplace transformation of both sides of Eq.(j22I) 



we get: 



vis) = ^ - - k[B]v{s)x{s) - ^v{s)Ms) - ^v{s)Ms) (23) 



(m — l)(g — 1 
function is given explicitly by 



where ~ = ^eg, which implies that the Laplace-transformed v{t) 



v(s) = ^ ^. (24) 

s (l + k[B]x{s) + ^ + ^^I^ + kv,,) 

Now, in virtue of equation ()12|1 . v{s) and the Laplace-transformed survival probability 
are related to each other as 



Pis) = (25) 
Hence, the Laplace-transformed survival probability obeys: 



Pis) 



S + 



. -1 



(26) 



l + k[B]xis) + ^Ms) + ^1^2is) 
It can be shown that xis), K,iis) and /?2(s) can be replaced by their value for s = 

to obtain the long time behaviour of Pit). Thus, inverting the latter equation in the limit 
t ^ oo, the A particle survival probability P(t) decays exponentially as: 
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m - exp ( . ,.^...:^"^1X;^MC1..(0) I (27) 



k[B][C]t 
k[B]m + ^ + -yr. 

and consequently, comparing the latter equation against the conventional form 



P{t)c^e^^{-Kff[B]t) (28) 

we have that the effective, overall reaction rate fee// describing the kinetics of diffusion- 
limited catalytically-activated reactions obeys: 

^ ^ ^ [B]m kM_ , ..q^ 
Kff k[C] [C] V'^[C] F'2 ^ 

which can be thought of as some "law of addition of inverse resistivities" and resembles 
(although has a more complex form) the classical result of Collins and Kimball for diffusion- 
limited trapping reactions with finite elementary reaction constant k js^. 

Equation is the central result of our analysis. Functions x(0),ft;i(0) and ^2(0) are 
studied in the next section. 



IV. RESIDENCE TIMES 

In this section we show that functions x(0),/ti(0) and /t2(0) have an apparent physical 
interpretation in terms of different residence times of Brownian paths in finite domains and 
may be evaluated in explicit form. 

A. One- particle's residence time 

The function £(0) entering the effective reaction rate is defined by: 

^(0) = V^l j drAdv\HAHAoG{vA, t; r° , 0). (30) 

One notices now that £(0) can be interpreted as the total time spent in a sphere of radius 
i? by a Brownian particle A, which started its diffusion at time t = at position (see 
figure El), averaged over all initial positions inside this sphere. In other words, this time is 
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the cumulative residence time inside the sphere up to an infinite observation time. In one 
or two dimensions it is infinite since the particle is certain to come back to the sphere, but 
in three dimensions it is finite since the particle can travel to infinity and thus is not certain 
to return. 

This residence time is well known 4l|, and can be calculated rather straightforwardly. 
Indeed, integrating first over the time with the change of variable y = -, we obtain: 

m = ^ (31) 



Now, since 



27r /-TT ^ 

TT-r sin 9 dOdcj) = tt--, (32) 

Jo ||rA-r^|| max{rA,r'^) 



we have: 



m 



1 An 



rdr^rl rdr\{r\f (33) 



Splitting next the second integral into the sum / = / + / , we get: 

Jo Jo Jta 

1 167ri?5 

which yields, eventually, the following result: 

£(0)-^ = hnDAR. (35) 

We hasten to remark that this expression, up to a numerical factor 5/6, coincides with 
the famous expression for the Smoluchowski reaction constant in three dimensions, Ks = 
AtiDaR. 

B. Two- particle's joint residence time 

Now, we turn to two other functions - Ki(0) and ^2(0), entering equation ()41|). The first 
one is formally defined by: 
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xHaoHboG{va, t- , 0)G'(rB„ t; 4, 0). (36) 

One may now notice that Ki(0) corresponds to the joint residence time in a sphere of radius 
R of particles A and -B, averaged over all initial positions of particles A and B inside the 
sphere (see figure E}. 

Such a residence time has been amply discussed in another work 42|, a nd here we will 
merely present the result of these calculations. It has been shown in Ref.j42] that /ti(0) is a 
complicated function of the diffusion coefficients Da and Db and is given explicitly by: 



^i(O) = — (— m (—\+—m(^\\ (37) 
' 207r \Da \Db Db \Da' ^ ^ ' 



where 



, , 2-101n(l + x) ln(l + x) , 7 

mix) = - 2 ^ ' + 16 arctan - -x^^^. (38) 

X ' x^^i 2 

Finally, the last undetermined function ^2(0) obeys: 



'^2(0) = yj^ J dt J d\idrAdY\drB,dr%^HAHB,HAoHB<i 

xG(r^,t;r° +u,0)G(rB.,t;r°,^ +u,0). (39) 

One notices that it may be interpreted as the joint residence time of particles A and B 
inside a sphere when the particles initially start from a sphere separated by the vector u 
from the residence sphere (see figure El). Note that the joint residence time ^2(0) is obtained 
when summing over all u. 

It is shown in Appendix that: 

T/'2 00 c; 

= 47rR{DA + Db)— ~ 1.15 x A7tR{Da + Db). (40) 
K2(,U) o34 
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C. Effective reaction rate 



Eventually, summing up the results of this section, we present an explicit expression for 
the effective reaction rate describing the kinetics of diffusion-limited catalytically-activated 
trapping reactions in terms of a suitably extended Wilemski and Fixman approach: 



J_^J_ [B] 
keff ~ k[C] ^ IAttDaRIC] ^ 

^g{DA,DB)7i^R'[C] ^ W^^nRiDA + D^) ^^^^ 

where g^^{DA, Db) = |l57"^ (^) ~^ (^)}' ^^''^^ being defined in equation 
Note that the correlations that are neglected in the Wilemski-Fixman approximation 
cannot be estimated precisely, so that it is difficult to find how they affect the different 
terms included into the reaction constant. However, these terms should not be considered 
as successive corrections to the usual expression for a bimolecular reaction. In fact, the joint 
residence time of two molecules on a catalytic site cannot be treated as a perturbation of the 
residence time of one molecule on this site, except in some limit cases. Clearly, our results 
should be compared to experimental or numerical results in order to be discussed. 

Consider now the behavior of the effective reaction constant in Eq. (j¥T|) in several limiting 
cases. In non-catalytic systems, in which the A and B particles may react at any point, 
which corresponds to an evident situation with [C] oo (but k[C] = K is kept finite), we 
find from Eq.(jiH) that 



keff^ ^^^A^{Da + Db)R' ^ ^ 

which represents, up to a numerical factor 385/334 ^ 1.15 in the second term, the classical 
Collins and Kimball result describing the effective reaction rate for trapping reactions 
involving diffusive A and B particles in non-catalytic systems. Note that the numerical 
factor 385/334 comes from the description of the reactive process which differs between the 
Collins-Kimball approach and the present case. 

Next, we turn to a different trivial situation when the particles B are present in a great 
excess, i.e. [B] oo, such that their diffusion becomes irrelevant. In this limiting case we 
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find from Eq.()4ip tfie following result: 

5 47rZ?^i?[C] 



-eff 



6 [B] 

such that the A particle survival probability follows 



(43) 



Pit) - exp (^-^AttDaRIC]?! (44) 

This is, again, a standard Smoluchowski-type (up to a numerical factor 5/6) prediction for 
trapping A + C ^ C reactions with immobile traps C. 

Finally, some simple analysis shows that for sufficiently small Db the third term on the 
right-hand-side of Eq. (j4ip dominates and the effective reaction constant follows 

Kff^giDA,DB)'iT''R'[C] (45) 

Surprising feature of this result is that k^ff is proportional not to the first power of the 
reaction radius R, but to the fourth power of it! Curiously enough, this prediction coincides 
with earlier results obtained for trimolecular reaction of the form A + A + C P + C using 
a„ extended Co.Uns-Ki.bal. approach Q. TK. anomalous dependence Kas been eonfi..n.ed 
by Molecular Dynamics simulations in Ref.j^. 



V. FLUCTUATION-INDUCED LONG TIME BEHAVIOR 

As shown in the introduction for trapping reactions involving diffusive particles, some 
fluctuation states can change signiflcantly the kinetics of diffusion-limited, catalytically- 
activated trapping reactions. The particular systems with random placement of the catalytic 
subvolumens the long-time kinetic behavior is described by a stretched-exponential function, 
so that the usual kinetic laws do not hold. 

Suppose that the A particle is initially at the origin, the traps B are also initially uni- 
formly spread in the reaction bath and the immobile catalytic subvolumens C are randomly 
distributed in the reaction bath such that the closest to the origin (i.e. to the A particle) 
subvolumen is at the distance p apart from it. Then, the A particle survival probability is 
evidently bounded from below by 

Pit)>PiacXPAip,t), (46) 
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where Piac is the probabihty of having a spherical region of radius p completely devoid of 
the catalytic suhvolumens C and PA{p,t) is, again, the probability that the A particle will 
not leave this region up to time t; in these conditions, the diffusive A particle may meet 
the diffusive traps B many times but the reaction can not take place since the necessary 
ingredient of the elementary act - the presence of the catalytic subvolumen C, will not be 
fulfilled. 

Consequently, the A particle survival probability obeys 

Pit) > exp X exp (-^i) (47) 

This bound is valid for any p and we have to choose such value of it, which provides the 
maximal lower bound. Maximizing the right- hand- side of the latter equation, we find that 
the optimal p is p ~ {DAt/27r[C]Y^^ , which yields 

P{t) > exp (^-^(27r)2/5[C]2/5(Z}^t)3/5^ , (48) 

i.e. the law similar to the one in Eq.Q describing the long-time evolution of the survival 
probability of a particle diffusing in presence of immobile, randomly placed traps. 

Since the right-hand-side of Eq. (j48p decays slower than exponentially, we may infer that 
at long-times the kinetics of diffusion-limited catalytically-activated trapping reactions will 
be supported by such fluctuation states and will be described by a stretched-exponential 
function of time. The comparison of|3Hlwith the classical kinetic law given bv 1^ shows that 
in all conditions the bound EH should only be considered for exceedingly long times, so that 
the classical law usually holds. 



VI. CONCLUSION 



To conclude, in this paper we have studied the kinetics of the catalytically-activated 
diffusion-limited reactions in Eq.(PQ) in the special case when B particles remain unaltered 
by reactions, i.e. the case of diffusion-limited catalytically-activated trapping reactions. In 
order to obtain an effective reaction rate for such bi-catalytic reactions, we have developed 
an analytical approach based on the work of Wilemski and Fixman 3^, which allowed us to 
calculate analytically the effective reaction constant. We have shown that this effective reac- 
tion constant comprises several terms which may be interpreted in terms of the residence and 
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joint residence times of Brownian paths in finite domains. We have demonstrated that the 
effective reaction constant exhibits a non-trivial dependence on the reaction radii, the mean 
density of catalytic subvolumens and particles' diffusion coefficients. Finally, we have dis- 
cussed the impact of several fluctuation states giving rise to anomalous fluctuation-induced 
contributions to the long-time kinetic behavior in such systems. Except in this asymptotic 
case, however the usual kinetic laws hold with the effective reaction constant calculated pre- 
viously. These results can be very useful not only in the theory of heterogeneous catalysis, 
but also in biology, when reactions can only take place on specific sites. 
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VIII. APPENDIX: JOINT RESIDENCE TIME IN A DISTANT SPHERE 

In order to evaluate explicitly the function given inEHl we first integrate the two propa- 
gators over u: 



duG{rA + u, t\r% 0)G{rB + u, t|r^, 0) 




0)GD^+Dg{rB-rA,t\rl-r%0) 
= Gda+Db{^b - rA,t|r^ - Ta^O), 



(49) 



where 



X 



Db 



{ya + u) + 



Da 



(rs + u) 



X" 



Da + Db 
Db 



(r^ + u) + 



Da + Db 
Da 




(50) 



Da + Db 



Da + Db 



and is the Gaussian propagator associated to the diffusion coefficient D. 
Then, the integral defined in ()39|) attains the following form: 
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\/'%(0) 



Hp^,aoHp^,B^Gda+Db{Pb - Pa5^|0, 0) 



(51) 



where p^ = — and = r^ — r^. The calculation of this equation requires the analysis 
of integrals of the following type: 



J dr^^HAoHp^^AO 



which represent the overlapping area between two spheres (see figure IH)). 
This area can be straightforwardly obtained: 



A = 2 r drTiiR' -r') = 1{R- pa/2)\4R + pa). 

Jpa/2 ^ 



Next, we have to calculate the integral 



(52) 



(53) 



\/'%(0) = / dtdpAdpsH{2R - \Pa\)H{2R - \p^\) 
Jo 

^{R- Pa/2)\4R + Pa)] [^iR - Pb^^^R + Pb) 



(54) 



Using the explicit expression for the propagators and integrating them over the time variable, 
we obtain 



and 



/ = / dpAdpsHi2R~\pA\)Hi2R-\p^\)- 



n m 

PaPb 



\Pb - Pa\ 



(55) 



2-77 /'TT 



Jo WPb ~ Pa\ 



1 An 

sin 9Bd6Bd(j)B 



max{pA,pB) 



(56) 



which yields 
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io io max{pA,pB) 

Consequently, the latter equation enables us to evaluate an explicit expression for K2(0): 



A-kR{Da + Db)— ~ 1.15 X A-kR{Da + Db). (58) 



^2(0) ' '334 
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FIG. 6 
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Figure legend 

Figure 1. Schematic representation of tlie bi-catalytic reaction: volume V, comprising 
a single A particle, m diffusing B particles and q immobile subvolumen C. 

Figure 2. Trimolecular reaction: the trajectories in bold type are those for which the 
reaction takes place, i.e. when one particle A and one B are jointly present in a catalytic 
domain. 

Figure 3. Residence time of a Brownian particle in a sphere. 

Figure 4. Joint residence time in a sphere. Blod lines: trajectories during common 
residence time. 

Figure 5. Joint residence time in a distant sphere. 

Figure 6. Overlapping area between two spheres. 
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